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Abstract 

The first-order differential equation of exponential relaxation can be gen- 
eralized by using either the fractional, derivative in the Riemann-Liouville 
(R-L) sense and in the Caputo (C) sense, both of a single order less than I. 
The two forms turn out to be equivalent. When, however we use fractional 
derivatives of distributed order (between zero and 1), the equivalence is 
lost, in particular on the asymptotic behaviour of the fundamental solu- 
tion at small and large times. We give an outline of the theory providing 
the general form of the solution in terms of an integral of Laplace type 
over a positive measure depending on the order-distribution. We consider 
with some detail two cases of fractional relaxation of distributed order: 
the double-order and the uniformly distributed order discussing the dif- 
ferences between the R-L and C approaches. For all the cases considered 
we exhibit plots of the solutions for moderate and large times. 
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1 Introduction 



The purpose of this paper is to study two types of fractional generalization of 
the classical relaxation equation. One type uses the fractional derivative in the 
sense of Riemann and Liouville, the other in the sense of Caputo. In its uses we 
distinguish between single and distributed orders of fractional derivatives. 

The plan of the paper is as follows. In Section 2, we recall the relevant 
properties of the fractional relaxation equations of a single order /3 = (5\ e (0, 1], 
in which the fractional derivative is intended in the Ricmann-Liouvillc (R-L) 
sense and in the Caputo (C) sense. The two forms are shown to be equivalent 
and the common solutions corresponding to a few orders are plotted. In Section 
3 we consider two general cases of fractional relaxation where in the R-L setting 
the order 1 — /3, in the C-setting the order [5 is distributed according to a non- 
negative weight function p{[5). For the fundamental solutions of these equations 
we provide a general formula obtained by the Titchmarsh theorem on Laplace 
inversion. In virtue of this, these solutions appear as real Laplace transform of 
a positive spectral function, and hence they are completely monotonic functions 
for t > in analogy with the fundamental solution of the fractional relaxation 
equation of a single order. In Section 4 we consider two typical cases of weight 
function: the case of two distinct orders < (5\ < (i 2 < 1 and the case of 
uniform distribution of orders between zero and 1. For these cases, by using 
Tauberian theory we provide the asymptotic expressions of the fundamental 
solution near zero and near infinity, that show the different role played by the 
order-distribution in the R-L and C approaches. Finally, concluding remarks are 
given in Section 5. For the reader's convenience we briefly recall in Appendix 
the essentials of the Fractional Calculus useful for understanding the notions of 
fractional derivative in the R-L sense and in the C sense. 



2 Fractional relaxation of single order 

The classical phenomenon of relaxation in its simplest form is known to be 
governed by a linear ordinary differential equation of order one, possibly non- 
homogeneous, that hereafter we recall with the corresponding solution. Denot- 
ing by t > the time variable, u = u(t) the field variable, and by tD 1 the 
first-order time derivative, the relaxation differential equation (of homogeneous 
type) reads 

t £>V*) = -Au(t), i>0, (2.1) 

where A is a positive constant denoting the inverse of some characteristic time. 
The solution of (2.1), under the initial condition it(0 + ) = 1 , is called the fun- 
damental solution and reads 

u(t) = e~ Xt , t>0. (2.2) 

From the view-point of Fractional Calculus (for a short review see the Appendix) 
there appear in the literature two ways of generalizing the equation (2.1), one 
way using the R-L, the other using the C fractional derivative. Adopting the 
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notation of the Appendix for the two derivatives, see (A.5)-(A.6), and denoting 
by P\ the common fractional order, the two forms read respectively for t > 

t D l u{t) = -X tD 1 ' u{t) , 0<j8<1, (2.3) 

and 

t D% u(t) = -\u(t) , 0</3<l, (2.4) 

where now the positive constant A has dimensions [t]~^ . If we assume the same 
initial condition, e.g. u(0 + ) = 1, it is not difficult to show the equivalence of the 
two forms by playing with the operators of standard and fractional integration 
and differentiatiorfj]. 

By applying in Eqs. (2.3)-(2.4) the technique of the Laplace transforms for 
fractional derivatives of C and R-L type, see (A.13)-(A.15), we get the same 
result for the fundamental solution, namely 

5(s) = ^Ta ' (2 ' 5) 

that, with the Mittag-Leffler functional Ep, yields in the time domain 

u(t) = E {~Xt p ) , 0<j9<1. (2.6) 

We agree to refer to the equation (2.3) or (2.4) as the simple fractional relaxation 
equation in the R-L or C sense, respectively. 



1 Both Eqs (2.3)-(2.4) are equivalent to the Volterra integral equation (of fractional type) 

u{t) = u(0+)-X t J t3 u(t). 

For example, we derive the R-L equation (2.3) from the fractional integral equation simply 
differentiating both sides of the latter, whereas we derive the fractional integral equation from 
the C equation (2.4) by fractional integration of order p. In fact, in view of the semigroup 
property (A. 2) of the fractional integral, we note that 

tJ 13 t D%u(t)= tJ P tJ 1 -* 3 tD 1 u(t)= t J 1 tD 1 u(t) = u(t)~u(0+). 

In the limit P = 1 we recover the relaxation equation (2.1) with the solution (2.2). The reader 
interested to have more details on the two forms of fractional relaxation may consult, for the 
R-L approach, the papers by Hilfer, Metzler and Nonnenmacher, see e.g. 11011241 . whereas for 
the C approach the papers by Caputo, Gorenflo and Mainardi, see e.g. |6l 1131 [20) . 

2 Let us recall that the Mittag-Leffler function Ep(z) (P > 0) is an entire transcendental 
function of order 1//3 , defined in the complex plane by the power series 

°° k 

£ 3 ( z ):=y^_J _ /3>0, zeC. 

PK ' r (j3k + 1) 

fe=0 

For details on it we refer to [TT][l7][T3]|25][26]. We remark that for t > the function Ep^Xt' 3 ) 
preserves the complete monotonicity of the exponential exp(— Xt): indeed it is represented in 
terms of a real Laplace transform of a non-negative function, 

„ N 1 f co c- Tt Xr' 3 sin(/37r) 

Eg(-Xt^) = — l ^—4 — dr, t>0, 0</3<l. 

^ ' 7T J r A 2 + 2A r$ cos(/37r) + r 2 ' 3 ' - ' 

However it decreases at t — * oo like a power with exponent -/3: Ep(-\t 9 ) ~ t'l 3 /[XF(1 - j3)\. 

If p = 1/2 we have for t > 0: E 1/2 (-^Vi) = e erfc(A v / t) ~ l/(A v / ^i) as t -» oo, where 
erfc denotes the complementary error function, see e.g. pQ. 
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Figure 1: Fundamental solutions of the fractional relaxation of a single order 
Pi = 1/4, 1/2,3/4, 1. Top: linear scales; Bottom: logarithmic scales. 

In Fig. 1 we show the solution (2.6) for a few values of the order j3 = Pi, 
Pi = 1/4, 1/2, 3/4, 1, by assuming A = 1: in the top plate for the time interval 
[0,10] (linear scales), and in the bottom plate for the time interval [10 1 ,10 7 ] 
(logarithmic scales). In the bottom plate we have added in dotted lines the 
asymptotic values for t — ► oo in order to better visualize the power-law decay 
expressed by i _/3l /r(l — Pi) for the cases < Pi < 1, whereas the case Pi = 1 is 
not visible in view of the faster exponential decay. In both plates we have shown 
in dashed line the singular solution for the limiting case Pi = 0, stretching the 
definition of the Mittag-Leffler function to Eq(z) = 1/(1 — z), the geometric 
series, 




(2.7) 



3 Fractional relaxation of distributed order 



3.1 The two forms for fractional relaxation 

The simple fractional relaxation equations (2.3)-(2.4) can be generalized by 
using the notion fractional derivative of distributed orde^. We thus consider the 
so-called distributed order fractional relaxation equation or fractional relaxation 
equation of distributed order, in the two alternative forms involving the R-L and 
the C derivatives, that we write respectively as 

t DVi) = —X f P (f3) t D 1 -^u(t)dp, (3.1) 
Jo 

and 

f p(f3)tD^(t)df3 = -Au«,(i), (3.2) 
Jo 

subjected to the initial condition u(0 + ) = u*(0 + ) = 1, where 

p(f3) > , and ( p(f3) d[3 = c> . (3.3) 
Jo 

The positive constant c can be taken as 1 if we want the integral to be normal- 
ized. For the weight function p([3) we conveniently require that its primitive 
P((3) = J p((3')df3' vanishes at (3 = and is there continuous from the right, 
attains the value c at f3 = 1 and has at most finitely many (upwards) jump 
points in the half-open interval < f3 < 1, these jump points allowing delta 
contributions to p(/3) (particularly relevant for discrete distributions of orders) . 

Since for distributed order the solution depends on the selected approach 
(as we shall show hereafter), we now distinguish the fractional equations (3.1) 
and (3.2) and their fundamental solutions by decorating in the Caputo case the 
variable u(t) with subscript *. 

As in [T2] the present analysis is based on the application of the Laplace 
transformation with particular attention to some special cases. Here, for these 
cases, we shall provide plots of the corresponding solutions. 

3.2 The integral formula for the fundamental solutions 

Let us now apply the Laplace transform to Eqs. (3.1)-(3.2) by using the rules 
(A. 15) and (A. 13) appropriate to the R-L and C derivatives, respectively. In- 
troducing the relevant functions 

A(s) = s f p{[3)s- p d[3 1 (3.4) 
Jo 

3 We find a former idea of fractional derivative of distributed order in time in the 1969 book 
by Caputo [p], that was later developed by Caputo himself, 4 5 and by Bagley & Torvik, see 
2 . A basic framework for the numerical solution of distributed-order differential equations 
has been recently introduced by Diethelm & Ford 9 , Diethelm & Luchko 1 101 and by Hartley 
& Lorenzo [ [T51IT9] . 
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and 

B(s)= f p{P)s dp, (3.5) 
Jo 

we then get for the R-L and C cases, after simple manipulation, the Laplace 
transforms of the corresponding fundamental solutions: 

and 

=•<«> = A^fr 

We easily note that in the particular case p{(3) = 8(f3 — fix) we have in (3.4): 
A(s) = s 1 - 131 , and in (3.5): B(s) = s 131 . Then, Eqs. (3.6) and (3.7) provide the 
same result (2.5) of the simple fractional relaxation. 

By inverting the Laplace transforms in (3.6) and (3.7) we obtain the funda- 
mental solutions for the R-L and C fractional relaxation of distributed order. 

Let us start with the R-L derivatives. We get (in virtue of the Titchmarsh 
theorem on Laplace inversion) the representation 

1 f°° 

u(t) = / e~ rt lm{u{re™)} dr , (3.8) 

77 Jo 

that requires the expression of — Im{l/[,s + AA(s)]} along the ray s = re™ with 
r > (the branch cut of the function s _/3 ). We write 



A [rc t7V j = p cos(7r 7 ) + ipsin(nj) , (3.9) 

where 

p = p{r) = \A(re™)\, 

7 = 7 (r) = larg[A(re-)]. (3 ' 10) 
Then, after simple calculations, we get 

/•CO 

u(t) = / c H(r; A) dr , (3.11) 
Jo 

with 

H(r;X)= 1 - a * P ™^> , ><>■ (3^) 
7r r z — 2\r p cos(7T7) + \ z p z 

Similarly for the C derivatives we obtain 

1 f°° 

u„(t) = / e- r *Im{u« (rc 47r )} dr , (3.13) 

n Jo 

that requires the expression of — Im {B(s)/[s(X + B(s))]} along the ray s = re™ 
with r > (the branch cut of the function s 13 ). We write 

B (rc in ^j = p* cos(7T7») + ip. t sin(-7r 7 *) , (3-14) 
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where 



!p* = p*{r) 
7* = l*{r) 
After simple calculations we get 



TT 



B(re™)\, 
- arg \B (r c 



(3.15) 




c~ n K(r;X)dr 



(3.16) 



with 



K(r;X) = — 



1 Ap* sin(7T7 >t ) 




(3.17) 



7r r X 2 + 2X p* cos(7T7*) + 



We note from (3.11) and (3.16) that, since H(r; X) and K(r; X) are non-negative 
functions of r for any A € R + , the fundamental solutions u(t) and u*(t) keep 
the relevant property to be completely monotone. 

The integral expressions (3.11) and (3.16) provide a sort of spectral repre- 
sentation of the fundamental solutions that will be used to numerically evaluate 
these solutions in some examples considered as interesting cases. 

Furthermore, it is quite instructive to compute for the fundamental solutions 
their asymptotic expressions for t — > and t — > oo because they provide their 
analytical (even if approximated) representations for sufficiently short and long 
time respectively, and useful checks for the numerical evaluation in the above 
time ranges. 

To derive these asymptotic representations we shall apply the Tauberian the- 
ory of Laplace transforms. According to this theory the asymptotic behaviour 
of a function f(t) near t — oo and t — is (formally) obtained from the asymp- 
totic behaviour of its Laplace transform f(s) for s —> + and for s —> +oo, 



respectively. For this purpose we note the asymptotic representations, 
from (3.6): 





being A(s)/s « 1/A 



(3.18) 



and from (3.7): 



f 1 B(s) 



s — > 0+ , being B(s)<<X, 



«.(*) ~ \ J r 3 ' A -I 



s— >+co, being B(s)»X. 



(3.19) 



B(s) 
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4 Examples 



Since finding explicit solution formulas is not possible for the relaxation equa- 
tions (3.1) and (3.2) we shall concentrate our interest to some typical choices for 
the weight function p(/3) in (3.3) that characterizes the order distribution. For 
these choices we present the numerical evaluation of the Titchmarsh integral 
formula, sec Eqs (3.8)-(3.12) for u(t) (the R-L case), and Eqs. (3.13)-(3.17) for 
u*(t) (the C case). The numerical results are checked by verifying the match- 
ing with the asymptotic expressions of u(t) and as t — > and t — > +oo, 
obtained via the Tauberian theory for Laplace transforms., according to Eqs. 
(3.18)-(3.19). 



4.1 The double-order fractional relaxation 

We now consider the choice 

p(P) =1*6(0-13!) +1*5(0-132), 0</3i</3 2 <1, (4.1) 

where the constants pi and P2 are both positive, conveniently restricted to the 
normalization condition pi + p 2 ~ 1. Then for the R-L case we have 

A(s)=p 1 s 1 - 01 +P2S 1 - 02 , (4.2) 

so that, inserting (4.2) in (3.6), 

^ (S) = S [l + \(p lS -^+p 2S -^)] ' (43) 
Similarly, for the C case we have 

B(s) = Pl s 01 +P2S 02 , (4.3) 

so that, inserting (4.3) in (3.7), 

~ I \ Pi S Pl + p 2 S 132 

uJs) = —r- 3 7TT ■ (4-4) 

V ^ S[\+P!sfc +P2sfo] V ' 

We leave as an exercise the derivation of the spectral functions H(r; A) and 
K(r; A) of the corresponding fundamental solutions, that are used for the nu- 
merical computation. The numerical results arc checked by their matching with 
the asymptotic expressions that we evaluate by invoking the Tauberian theory 
and using Eqs. (3.18)-(3.19) jointly with Eqs (4.2)-(4.3) respectively. 

For the R-L-case we note that in (4.2) s 1- ^ 1 is negligibly small in comparison 
with s 1_ ^ 2 for s — ► + and, viceversa, s 1- ^ 2 is negligibly small in comparison to 
s 1 -0i f or s _ > +oo. Similarly for the C-case we note that in (4.3) s 132 is negligibly 
small in comparison to s^Tor s — > + and, viceversa, s^ 1 is negligibly small in 
comparison s^for s — ► +oo. 
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As a consequence of these considerations we get for the R-L case, if fa < 1, 

+ , 

u(s) ~ { "P* (4.5) 




so that 



+00 , 



£ — ► +00 . 



„«,. < A '* m (4.6) 

1 - Api — , t -> 0+ . 





Figure 2: Fundamental solutions of the R-L- fractional relaxation of double or- 
der in some {fix, fa) combinations: {1/8,1/4}; {1/4,1/2}; {1/2,3/4}; {3/4,1}. 
Top: linear scales; Bottom: logarithmic scales. 



We note that Eq. (4.5a) and henceforth Eq. (4.6a) lose their meaning for 
fa = 1. In this case we need a more careful reasoning: we consider the expression 
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for s — > provided by (3.18) as it stands, that is 
In virtue of the Laplace transform pair 

see Eq. (1.80) in [25] . where Ea tV denotes the Mittag-Leffler function in two 
parameter^ we get, with q — P2/(Api) aud ji = v = 1 — /3±, as t — > +00 : 

«(*) - t" - *"* ^- A>1 -ft (-gt 1 -*) = -t 1 - ^i-ft • (4-9) 

Api Api at 

Taking into account the asymptotic behaviour of the Mittag-Leffler function, 
we finally get 

<t )~\^±=£t-( 2 -M as t^+00. (4.10) 



Similarly for the C case we get: 



~ { 1 /. A (4-11) 



so that 



Pi 




A 






'i-A. 










pi 





H 



a r(i-/3i) ' 



i — » +00 . 



" I'l - < i ^ (4-12) 

r(i + /3 2 ) ' ^ 0+ - 

We exhibit in Figs. 2-3 the plots of the fundamental solutions for R-L and C 
fractional relaxation, respectively in some {^1,^2} combinations: {1/8,1/4}; 
{1/4,1/2}; {1/2,3/4}; {3/4,1}. We have chosen pi = p 2 = 1/2 and, as usual 
A = 1. From the plots the reader is expected to verify the role played by 
the different orders for small and large times according to the corresponding 
asymptotic expressions, see Eqs. (4.6), (4.9)-(4.10) and (4.12). 



4 The Mittag-Leffler function E^.y^z) (5R{/i} > 0, v £ C) is denned by the power series 

°° k 

It generalizes the classical Mittag-Leffler function to which it reduces for v = 1. It is an 
entire transcendental function of order l/5R{/x} on which the reader can inform himself by 
again consulting [HI El [ED 125] [26] . With fj,, v 6 R the function E^^(-x) (x > 0) turns a 
completely monotonic function of x if < fj, < 1 and v > fi, see e.g. 1231 . This property is 
still valid when x = qt^ (q > 0). In particular, for < fi = u < 1 we note 

qt^ 1 -^ (-qn = -±E» (-qt») ~ — ft f^+D , t - +CO . 

at g I (1 — fj,) 
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Figure 3: Fundamental solutions of the C-fractional relaxation of double order in 
some {/3i,/3 2 } combinations: {1/8,1/4}; {1/4,1/2}; {1/2,3/4}; {3/4,1}. Top: 
linear scales; Bottom: logarithmic scales. 



4.2 The uniformly distributed order fractional relaxation 

We now consider the choice 



<P< 1. 



(4.13) 



For the R-L case we have 




(4.14) 



hence, inserting (4.14) in (3.6) 



logs 



(4.15) 



slogs + A (s — 1) ' 
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For the C case we have 



B(s) = [ s^dp 
Jo 



logs 



hence, inserting (4.16) in (3.7), 



1 



u*(s) = - 



s- 1 



s A log s + s — 1 



1 1 A logs 
s s A log s + s — 1 



(4.16) 



(4.17) 



We note that for this special order distribution we have A(s) — B(s) but the 
corresponding fundamental solutions are quite different, as we see from their 
Laplace transforms (4.15) and (4.17). 

Then, invoking the Tauberian theory for regularly varying functions (power 
functions multiplied by slowly varying function^, a topic adequately treated in 
the treatise on Probability by Feller [12], Chapter XIII. 5, we have the following 
asymptotic expressions for the R-L and C cases. 

For the R-L case we get 



u(s) 



logs 
A(s-l)' 



1-A- 



s log s 



0+, 
+oo . 



(4.18) 



so 



u(t) 



A 6 flW ~A7< « 



1 



A 



|log(l/t)| 



i^0+. 



(4.19) 



In (4.19a) £ x (t) := 

and the Laplace transform pair (29.3.100). 
For the C case we get 



■ du denotes the exponential integral, see [T], Ch. 5 



1 



u*(s) 



Aslog(l/s) 
1 A logs 



0+, 



(4.20) 



so 



«*(*) ~ { A \ogt ' 

1 - Ailog(l/t), i^0+. 



(4.21) 



5 Deflnition: We call a (measurable) positive function a(y), defined in a right neighbour- 
hood of zero, slowly varying at zero if a(cy) / a(y) — » 1 with y — » for every c > 0. We call a 
(measurable) positive function b(y), defined in a neighbourhood of infinity, slowly varying at 
infinity if b(cy) /b(y) — > 1 with y — » oo for every c > 0. Example: | logt/| 7 with 7 £ R . 
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10 1 1(f 10 3 10 4 10 5 10 6 10 7 

Figure 4: Fundamental solutions for R-L and C uniformly distributed fractional 
relaxation in comparison with some of them for single orders. Top: linear scales; 
Bottom: logarithmic scales. 

In Fig. 4 we display the plots of the fundamental solutions for R-L and C 
uniformly distributed fractional relaxation, adopting as previously, in the top 
plate, linear scales (0 < t < 10), and in the bottom plate, logarithmic scales 
(10 1 < t < 10 7 ). 

For comparison in the top plate the plots for single orders /3i = 0, 1/2, 1 are 
shown. We note that for 1 < t < 10 the R-L and C plots are close to that for 
(3\ = 1/2 from above and from below, respectively. 

In the bottom plate (where the plot for /?i is not visible because of its faster 
exponential decay) we have added in dotted lines the asymptotic solutions for 
large times. We recognize that the C plot is decaying much slower than any 
power law whereas the R-L plot is decaying as t~ x ] this means that for large 
times these plots are the border lines for all the plots corresponding to single 
order relaxation with /3i 6 (0, 1). 
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5 Conclusions 



We have investigated the relaxation equation with (discretely or continuously) 
distributed order of fractional derivatives both in the Riemann-Liouville and in 
the Caputo sense. Such equations can be seen as simple models of more general 
distributed order fractional evolution in a Banach space where the relaxation 
parameter A is replaced by an operator A acting in this space. A relevant ex- 
ample is time-fractional diffusion where in the linear case the individual modes 
exhibit fractional relaxation. Our interest is focused on structural properties of 
the solutions, in particular on asymptotic behaviour at small and large times. In 
both approaches we find that the smallest order of occurring fractional differen- 
tiation determines the behavior near infinity, but the largest order the behaviour 
near zero, in analogy to the special form of time-fractional diffusion explicitly 
governed by the distributed order derivative as in [8] and in [T51 [2TJ [27] . We 
see that the two parameters (3i and fa play opposite roles in our two cases (R-L) 
and (C). The topic deserves further study in several directions, e.g. in in terms 
of integral transforms and special functions like those of Mittag-Lefner type. 
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Appendix: Essentials of fractional calculus 

For a sufficiently well-behaved function f(t) (t € R + ) we may define the frac- 
tional derivative of order fj, (m — 1 < fi < m, to e IN), see e.g. [T51 1^5] . in two 
different senses, that we refer here as to Riemann-Liouville (R-L) derivative and 
Caputo (C) derivative, respectively. Both derivatives are related to the so-called 
Riemann-Liouville fractional integral of order a > defined as 

f/°7(*) := pTT / (t-rr^fi^dr, a>0. (Al) 

r M Jo 

We recall the convention t J° = I (Identity operator) and the semigroup prop- 
erty 

t J a t jP= t J P tJ a = tJ a+P , a,p>0. (A.2) 

Furthermore 

tJ a f = J^ 1 } r * 7+a , «>0, 7 >-l, t>0. (A3) 
T(7 + 1 + a) 

The fractional derivative of order /i > in the Riemann-Liouville sense is defined 
as the operator t_D M which is the left inverse of the Riemann-Liouville integral 



14 



of order /i (in analogy with the ordinary derivative), that is 



(A4) 



If m denotes the positive integer such that m — 1 < fi < m , we recognize from 
Eqs. (A.2) and (A.4) t D» f(t) := t D m t J m ~* f{t) , hence 



i 



r(m - n) J (t - t) 



/(r)dr 



1 < /i < m, 



(A5) 



/i = m. 



For completion we define t D° = / . On the other hand, the fractional derivative 
of order fx > in the Caputo sense is defined as the operator t Z?* such that 
tDH fit) := t J m -» t D m f{t) , hence 



1 



lD , l[t )=\ r(m- M ) 



* f (m) {r)dT 



1 < /i < m, 



(A6) 



:/(*)• 



/i = m. 



Thus, when the order is not integer the two fractional derivatives differ in that 
the derivative of order m does not generally commute with the fractional inte- 
gral. We point out that the Caputo fractional derivative satisfies the relevant 
property of being zero when applied to a constant, and, in general, to any power 
function of non-negative integer degree less than m , if its order \x is such that 
m — 1 < ji < m . Furthermore we note that 



r( 7 + i) 
r( 7 + 1 - m) 



f-^, ^>0, 7>-l, t>0. 



(A.7) 



Gorenflo and Mainardi [13] have shown the essential relationships between the 
two fractional derivatives (when both of them exist), 



k= 
m— ] 



fc=0 



k 
k\ 

u 

fW(0+)t k -v 
T(k-fi + l) 



1 < fi < m . 



(A8) 



In particular, if m — 1 we have 



t D» [f(t) - /(0+)] , 
Jy> T{l-n) 



< fi < 1 . 



(A9) 



The Caputo fractional derivative, practically ignored in the mathematical trea- 
tises, represents a sort of regularization in the time origin for the Riemann- 
Liouville fractional derivative. We note that for its existence all the limiting 
values / < - fc - l (0 + ) := lim fit) are required to be finite for k = 0, 1, 2 m — 1. 
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We observe the different behaviour of the two fractional derivatives at the 
end points of the interval (m — 1, m) namely when the order is any positive 
integer: whereas t D^ is, with respect to its order \i , an operator continuous at 
any positive integer, t D* is an operator left-continuous since 



lim t DU{t) = f( m - 1 \t)-f( m -V(Q+), 

/i— »(m— 1) + 

lim t DU(t) = f {m) (t)- 



(AlO) 



We also note for m — 1 < /z < m , 



t D» f(t) = t D» g(t) f(t) = g(t) + ]T Cj , (All) 

3=1 
m 

t D% f{t) = t D» g(t) ^=> f{t) = git) + 9 . (A12) 

3=1 

In these formulae the coefficients Cj are arbitrary constants. Last but not least, 
we point out the major utility of the Caputo fractional derivative in treating 
initial-value problems for physical and engineering applications where initial 
conditions are usually expressed in terms of integer-order derivatives. This can 
be easily seen using the Laplace transformation, according to which 

m — 1 

4 = ^/(5)- ^ S ^ 1 - fc / (fe) (0 + ), m-l<n<m, (A13) 



oo 



where f(s) = £{/(*); s} = / e " s< fit) dt , s G C, and / (fc) (0 + ) := lim /(*)• 

Jo . *^o+ 

The corresponding rule for the Riemann-Liouville derivative is more cumber- 
some: for m — 1 < /i < m it reads 



m—1 



£ /(t); s} = f(s) - [t * tJ {m -^] f(0 + ) s"- 1 -' , (A14) 



fe=0 



where, in analogy with (A. 13), the limit for £ — ^ + is understood to be taken 
after the operations of fractional integration and derivation. As soon as all the 
limiting values /^ fe ^(0 + ) are finite and m — 1 < (J, < m, the formula (A. 14) 
simplifies into 

£{ t D»f(t); S } = S »f(s). (A15) 

In the special case /(*)(0+) = for k = 0,1, to- 1, we recover the identity 
between the two fractional derivatives, consistently with Eq. (A. 8). 

For more details on the theory and applications of fractional calculus we 
recommend to consult in addition to the well-known books by Samko, Kilbas 
& Marichev [25], by Miller & Ross [55], by Podlubny those appeared in 
the last few years, by Kilbas, Srivastava & Trujillo [T7], by West, Bologna & 
Grigolini [55], and by Zaslavsky [55] . 
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